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The transmission of vector infectious diseases, which produces complex spatiotemporal patterns, 
is analyzed by a periodically forced two-dimensional cellular automata model. The system, which 
comprises three population levels, is introduced to describe complex features of the dynamics of the 
vector transmitted dengue epidemics, known to be very sensitive to seasonal variables. The three 
coupled levels represent the human, the adult and immature vector populations. The dynamics 
includes external seasonality forcing (rainfall intensity data), human and mosquito mobility, and 
vector control effects. The model parameters, even if bounded to well defined intervals obtained 
from reported data, can be selected to reproduce specific epidemic outbursts. In the current study, 
explicit results are obtained by comparison with actual data retrieved from the time-series of dengue 
epidemics in two cities in Brazil. The results show fluctuations that are not captured by mean-field 
models. It also reveals the qualitative behavior of the spatiotemporal patterns of the epidemics. In 
the extreme situation of absence of external periodic drive, the model predicts completely distinct 
long time evolution. The model is robust in the sense that it is able to reproduce the time series of 
dengue epidemics of different cities, provided the forcing term takes into account the local rainfall 
modulation. Finally, the dependence between epidemics threshold and vector control undergoes a 
transition from power law to stretched exponential behavior due to human mobility effect. 

PACS numbers: 87.18.-h, 87.16.aj, 87.19.xd, 87.15.A 



I. INTRODUCTION 

Understanding the rather complex dynamics of trans- 
missible diseases is of utmost importance for improving 
life quality, and even the survival of some human popula- 
tion groups. To achieve this, interdisciplinary efforts are 
necessary, which certainly include the use of the recently 
techniques developed to study complex systems [UEIE]. 
At the beginning of 21 century, both directly transmitted 
diseases, like tuberculosis and AIDS, as well as vector- 
transmitted diseases, such as dengue and malaria, are 
still not controlled. In modern life, the intense flux of 
people at global level and within large cities [I] increases 
the complexity of the propagation of transmitted diseases 
[5 . For vector-transmitted diseases, there are already 
indications that climatic conditions and vector mobil- 
ity may increase the number of cases [5J. In the case 
of dengue, an arboviral disease transmitted to humans 
by Aedes mosquitoes (mainly Aedes Aegypti), several de- 
terminant factors for its transmission are found in large 
urban centers [7]: human concentration, large inter- and 
intra-city human mobility, the climatic conditions for the 
vector proliferation (high humidity and temperature be- 
tween 15°C and 40°C). Accordingly, it is found that the 
dengue outbursts are quite sensitive to seasonal varia- 
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tions in pluviometric precipitations, humidity and tem- 
perature. The disease, which may be caused by four 
different virus serotype (DenVl-DenV4), reaches yearly 
some 50 millions people in more than 60 countries, with 
- 21000 casualties [8]. 

Since 1992 [3], ordinary differential equation (ODE) 
models have been proposed to analyze dengue inter-host 
dynamics and the effect of vector control actions. More 
recently, some attempts to introduce the spatial depen- 
dence on the disease propagation have been reported, 
using both partial differential equation (PDE) [TU] and 
cellular automata (CA).llJ models, and other data anal- 
ysis techniques [T2] . In [TT] , the authors proposes a model 
that takes into account only the description of mosquito 
population, which may be found in the adult phase, and 
the immature phase comprising several stages. However, 
a more accurate description of the dengue propagation 
must include, besides the interaction among these pop- 
ulation groups, the vector mobility, effect of control ac- 
tions, and an explicit climatic periodic forcing on the 
population variables. To our knowledge, no previous in- 
vestigation has taken into account all of these factors. 

In this work, we investigate an inter-host three level 
CA model, which describe the pertinent population 
groups in a urban environment: human, adult vector 
mosquito, and immature vector in the aquatic phase. As 
we will detail later on, it includes all of the quoted effects: 
external forcing to describe the environment influence on 
the vector life cycle, as well as other interaction terms de- 
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scribe the effect of human and vector mobility and control 
actions. The results provided by the model reproduce 
actual time series from some well document dengue epi- 
demics in specific years urban centers in Brazil. Besides 
that, they also qualitatively agree with main features of 
the spatiotemporal patterns. We also show that, in the 
absence of a periodic forcing, the actual epidemic out- 
bursts are not reproduced, supporting the claims of the 
importance of climatic aspects in the triggering of local 
events. Finally, as the model describes the behavior of 
the exposed population for larger time intervals under 
the presence of climatic seasonal variations, it is possible 
to follow the effect of vector control actions. In such case, 
our results indicate a power law dependence between the 
epidemic threshold and the parameter describing the in- 
tensity of vector control. 

The current description of vector-transmitted diseases 
goes along several successful works based on CA intra- 
host disease propagation models (for instance, AIDS, 
[13] , malaria [14] , cancer [15] ) and also on inter-host mod- 
els [TB]. It is also worth mentioning that the presence of 
multiple CA interacting levels in epidemic models has 
been explored in alternative topologies, as that of com- 
plex networks where nodes represent patches of regular 
lattices [T7J submitted to a contact process dynamics [IB] ■ 

The paper is organized as follows: in Section [IT] we 
introduce the CA local rules, comparing them to other 
models in the literature. Section |lll] discusses the choice 
of parameter values in our simulations. In Section |IV| 
we present our results, comparing them with actual data: 
the simulated time series (IV A) resulting from the peri- 



odic forcing seasonal effects, the simulated spatiotempo- 
ral patterns and the vector control associated to human 
mobility effect. Finally, Section [V] closes the paper with 
concluding remarks and perspectives. 



II. THE MODEL 

Some of the basic interaction mechanisms and external 
effects be included in our three level CA model have been 
used, in other context, by previous ODE models reported 
in the literature. The first attempt [9J considered a com- 
partment model, in which humans follow SEIR (suscep- 
tible, exposed, infected and removed) dynamics. Since 
mosquitoes usually die before being removed, the authors 
consider that they follow a simpler three-compartment 
SEI version. On the other hand, climatic effects were 
modeled by seasonal variations of model parameters by 
an ODE system |19j . Tuning models by comparison to 
actual data have also been attempted, e.g., by the esti- 
mation of the basal transmission rate for age-stratified 
data from Thailand [5U]. Other models have considered 
the role of a unique vector in the transmission of multiple 
diseases, as more than one dengue serotype [2TJ [23] [24] 
or the concurrent transmission of yellow fever in dengue 
infested areas [25 . Finally, the effect of vector control 
have already been explicitly analyzed in ODE models 
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Each of three CA levels consists of a two-dimensional 
square lattice with N s = L x L sites. Correspondingly, 
the CA is subjected to closed boundary conditions be- 
cause it mimics dengue transmission in a city. If we 
compare the results to actual data, each neighborhood 
corresponds to a set of distinct spatial units (census sec- 
tors) into which the reported cases are assigned to. Each 
site in the distinct levels describe, respectively, the local 
populations: human (H), mosquito (M) and immature 
vector in the aquatic phase (A). The CA inter-layer in- 
teraction rules couple, locally, the three involved levels 
due to the interactions between H and M levels, and the 
A to M flux of the vector population. The CA Moore 
neighborhood with radius 1 allows, for each site of a given 
layer, a maximum of 9 neighbors in the level it interacts 
with (see Figure [TJ). We restrict ourselves to the one- 
serotype situation, although the model can be extended 
to simulate the dynamics with more than one serotype. 




FIG. 1: Diagram of different lattices: humans (H), 
mosquitoes (M) and aquatic phase (A). Note that each el- 
ement of lattices H and A 'sees' up to nine neighbors of the 
lattice M (and vice- versa). 

According to previously indicated models, in the A 
phase, the vector is found in one of 4 compartments: 
egg (E), larvae (L), pupae (P) and breeding (B). The 
M phase comprises 3 compartments: susceptible (SM), 
exposed (EM), and infectious (IM). Finally, consider- 
ing only one serotype, there are 4 possible compartments 
for H sites: susceptible (SH), exposed (EH), infectious 
(IH), and recovered (RH). Moreover, sites of A and 
M levels can be in empty states, denoted by EAS and 
EMS. The local interaction rules, based on the ento- 
mological [27] and epidemiological aspects [SS], are such 
that, for each level: (see Figure [2]). 

A level: E, L, P and B states evolve from the pre- 
ceding one after the E eclosion period t e , L phase 
period t; and P phase period t p . An empty site 
EIM may be replaced with probability f s (t) by an 
E state, if there is at least an occupied site in its 
Moore neighborhood at the M level. The transi- 
tion from E to L compartments also depends on 
f s (t), much as the persistence of B, which releases 
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fi, which is the same for all sites, is to globally increase 
the SM -> EM and SH -»■ EH probability transitions, 
without any influence from the neighborhood population 
in the other level. 

Finally, the decrease of populations in M level result- 
ing from vector control actions is included by the follow- 
ing additional rule: the natural M death probability is 
increased by an additional amount p a d m , which reduces 
the adult mosquitoes on any state of the M level. 
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FIG. 2: A schematic representation of the local rules of the 
model. 



an adult mosquito SM to a EPM site of the M 
level. 



M level: The population in the M level results 
from the dynamics in A phase. Adult population 
M dies according to a death probability pdm in any 
state. The transition from a SM site into EM de- 
pends on the number of IH sites in its neighbor- 
hood in the H level, on the local effective biting 
humans- mosquitoes probability p e hm, and on the 
human mobility \x. An EM site becomes IM after 
the M virus latent period U m . 



H level: In a similar way to the SM — > EM tran- 
sition, a SH site changes to EH according to the 
local effective mosquitoes-humans biting probabil- 
ity (pemh), the number of IM sites in the M level 
neighborhood, and on the human mobility \i. EH 
becomes infectious IH after the H virus latent pe- 
riod tih, and IH becomes recovered RH after the 
viremia period t v . 

Note that, in the above level descriptions, we already 
included relevant features of dengue transmission that 
we have called the attention in Section 1. Seasonal in- 
formation (rainfall intensity) is used as input data [5] [7] 
by tuning the time dependence of the f s (t) probability, 
using a Fourier expansion of the actual rainfall series. 
If the time series do not include daily entries, or is not 
complete over the whole simulation period, interpolation 
or addition of random noise to the day average taken 
over a few years can be used. Global infection probabil- 
ities between H and M populations, due to mobility in 
private and public transport systems, is described by a 
global (mean-field) mobility parameter [i. The action of 



III. PARAMETER VALUES 

The CA parameters introduced in the previous section 
can be classified into four classes, according to the indi- 
vidual process they describe: 1) Spatial parameters, as L 
and /j,; 2) Temporal parameters: t e , tj, t p , ti m , t v , and tin] 
3) The probabilities of transmission and mosquito death 
parameters: p emh , p e hm, an d p md ; 4) Vector control pa- 
rameter: Padrn- 

The values of spatial parameters are obtained by tak- 
ing into account the data of a given urban center. We 
estimate the size L of the lattice (number of sites = I 2 ) 
by the area of the city (A c ) and the flight radius of the 
vector (R). More specifically, we assume that A c = L 2 a, 
where a is the area of one cell, while R corresponds to 
the average (Moore) neighborhood radius. This way, we 
have 



R = 



y/a(l + >/2) 
2 



L = 



^(1 + 
2R 



(1) 



As the dispersion of Aedes aegypti due to its flight 
rarely exceeds 100m [33], we assume R = 100m. The 
range of values of \i was estimated by requiring that the 
model reproduces the same behavior of the histogram of 
the number of census sectors with, at least, one reported 
dengue case during the corresponding time period. 

We assumed fixed values (within the range presented 
in Table [I]) for the probabilities of transmission p e h m — 
Pemh = 0.75 [H] and of mosquito death p m d = 1/7 = 
0.143. For vector control parameter, when is the case, 
we scrutinize the complete interval from to 1. 

Choosing the CA iteration time unit to be one day, 
we are able to set value intervals for several temporal 
parameters according to the literature (see Table [l} . To 
obtain baseline values for temporal parameters and epi- 
demic threshold, we adapt the epidemiological definition 
of an epidemic process [33] to our model simulations. A 
disease is considered epidemics if the annual incidence I, 
the number of reported case to susceptible population, is 
above a certain (epidemics) threshold I t h- Therefore, Ith 
may be given by 



Ith =<!> +2(J, 



(2) 



where the average incidence < I > is calculated with 
respect to the last N years and a corresponds to the 



4 



standard deviation. To obtain corresponding model val- 
ues, we run the program for N different random seeds. 
We recall that, as for actual cases of vector transmit- 
ted diseases, several numerical simulations resulting from 
different random seeds die out in the first weeks, being 
characterized as small endemic processes. 

After the evaluation of < 7 > and <r, we run the pro- 
gram as many times as necessary to get K independent 
samples with 7 > I t h- Although we perform the numeri- 
cal simulations of the model for large time intervals, our 
analysis can be restricted to 364 time unit intervals if we 
want to compare the results with actual data of one year 
epidemics series. The output data are the time series 
of density of each state in the 77, M and A levels of the 
CA model, and the spatiotemporal configurations at any 
time step. The cpu time increases according to 7 3 and 
linearly with the number of samples. 

Finally, based on the range of values in Table [T]for tem- 
poral parameters, simulations have been conducted for an 
initial set of parameter values. Then, we investigate the 
effect of changing one by one parameter, while holding all 
the others fixed. This way, we identify the baseline values 
that minimizes the error between the actual time series 
and the simulated time series. We perform several tests 
in order to check the robustness of the chosen initial set 
of parameter values. For a systematic analysis of param- 
eter values, we considered an average of M simulations 
samples, identifying the best output for the purpose of 
comparison with one actual epidemics time scries. This 
is achieved by the analysis of the minimum discrepancy 
between actual and simulated time series: 



A. The seasonal effects: actual and simulated time 

series 

The 1995 weekly rain intensity and reported num- 
ber of new dengue cases Id (incidence) are shown in 
Figure [3] where the data have been normalized by the 
largest input for the sake of comparing the tendency of 
the curves. As, in this case, temperature and humidity 
are quite stable, rainfall is the most important climatical 
factor for dengue propagation. Indeed, the Pearson cor- 
relation varies from 0.49 to 0.76 for, respectively, weekly 
and monthly sampled data. As it will be clear from the 
discussion of our results, such increase in the correlation 
in value is due to a roughly two week delay time between 
the two signals. When the series are clustered in large 
time windows, such effects become much smaller. The 
daily rainfall data was provided by the Brazilian govern- 
ment |36j . 




where T is the number of days, a* is the actual incidence 
and Si is the simulated incidence of day i. 

Once estimated the baseline of temporal parameter 
values, the analysis of minimal discrepancy is also ap- 
plied to select the best sample in comparison to actual 
data. 



IV. RESULTS 

In order to validate the model, we consider the data of 
the first dengue epidemics (DenV-2) in 1995, Salvador, 
Brazil |35j . when its population p c = 2.3 million habi- 
tants distributed over an A c = 313 x 10 6 m 2 area. In 
1995, the average daily temperature was 25.89°C with 
1.47 standard deviation. The city yearly mean precipita- 
tion is 1980 mm/year, while seasonal effects concentrate 
precipitation in the months March- August. 



FIG. 3: Normalized time series of weekly Id in Salvador 
(black-squares) and normalized time series of weekly Tr 
(grey-circles) for 1995. The normalization factor are 846 cases 
and 373 mm respectively. 

The simulations are based on the function f s (t) corre- 
sponds to the Fourier expansion 

12 

f s (t) = a Q + ^2a j cos(Trjt/26) + b :j sin(irjt/26) (4) 

with a = 0.13585, a x = -0.12872, b x = 0.05071, a 2 = 
0.0502, b 2 = -0.0882, a 12 = 0.00744, b 12 = 0.04713. The 
total contribution of the the remaining coefficients (Zj, 6j 
can be neglected. 

The 10, 831 reported dengue cases in Salvador dur- 
ing 1995 were geo-referenced by epidemiological week 
(52 temporal units) and census sectors (2600 spa- 
tial units) [35 . Note that, due to large official sub- 
notifications (26 %), the actual number of cases is much 
larger. In Salvador, the epidemics peak occurs before 
the rainfall peak, what can be justified by the fact that, 
due to the high intense pluviometric precipitation peak, 
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the rainfall washes out the vector in the immature phase. 
As well will see later, this may not happen in other ur- 
ban centers. According to expression (JTJ, we are lead 
to the value L = 214. On the other hand, the value 
ji = 5 x 10~ 4 , has been selected from the interval where 
the model is able to reproduce the exponential behavior 
in the probability distribution of observed new cases in a 
year among 2600 sensus sectors (not shown). 

Assuming that there is one infected individual in each 
site of H lattice, the best sample is able to reproduce 
the actual data quite well, as shown by the normalized 
actual and simulated incidence time series in Figure [4] 
We normalize both the actual and simulated time series 
for the purpose of avoiding distortions due to large sub- 
notifications. To set up the correspondence between In 
and the simulated incidence IHjy, that is, the number 
of new infected humans at a time step, we use the scale 
factor L 2 /p c . The normalization factors for the actual 
and the best simulated incidence time series result, re- 
spectively 17 and 11. 




W«aks 



FIG. 5: Normalized average over 20 samples of simulated 
weekly incidence (IHn) when f s (t) = 1 (black square) and 
f s (t) = sin(2nt/52) (grey circle). The normalization factors 
for f s (t) = 1 and f B (t) = sin(2irt/52) are 35.5 and 17.3 re- 
spectively. The data has been smoothed by averaging on three 
consecutive weeks. Parameter values are the same as in Fig- 
ure m 




FIG. 4: Actual (Id) and simulated (IHn) weekly incidence 
time series of Salvador in 1995 normalized by largest single 
input. The data has been smoothed by averaging on three 
consecutive weeks. Circles, squares and triangles indicate, 
respectively, Id, the best individual sample, and average value 
over 20 samples taken from random seeds. The normalization 
factors for Id, the best IHn and the averaged IHn, are 17, 
11 and 14 respectively. Consider one sample and the following 
parameter values: t e = 5, ti = 5, t p = 3 , U m — 7, tin = 6, 
t v = 6, p ehm = p emh = 0.75, p md = 0.143, >i = 5x 10" 4 . 

Note that the delay between the peaks of Id and the 
best individual sample is much smaller that the delay 
between Tr and Id in Figure [3] even considering the 
averaging on three consecutive weeks which amplify the 
delay effect. Although this effect is also amplified for the 
average over some samples, it is still smaller than the 
delay between Id and Tr in Figure [3] 

To emphasize the importance of the periodic forcing to 
recover the reported Id values, we draw, in Figure|5] the 
time evolution according to two hypothetical scenarios. 



They were obtained by replacing f s (t) , in first place, by a 
constant value, and afterwards by a simple periodic sine 
function. The resulting incidence counts differ substan- 
tially from the typical patterns in Figure [5] The impor- 
tance of such external drive, which is a crucial aspect of 
vector transmitted diseases, has been neglected in most 
of analyzed models with time and space dependence. 
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FIG. 6: The normalized rain intensity (Tr), actual (Id) and 
simulated (IHn) incidence (by month) in Mossoro in 1999. 
The normalization factors are 71 cases, 149.5 mm, and 72 
cases respectively. Consider the best sample and the following 
parameter values: t e = 4, ti = 7, t p = 7 , U m — 6, tih = 5, 

tv = 6, Pehm = Pemh = 0.75, p m d = 0.143, fl = 1.0 X 10" 6 . 

The used lines-symbols are, respectively, dashed white-circle, 
solid black-square, and dotted grey-triangle. 

The importance of seasonal aspects for the observed 
dynamics can be further exemplified by running the 
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model with the data of other urban centers. For instance, 
we consider the 1999 dengue epidemics in Mossoro, in 
Northeast Brazil [35], for which rainfall peak precedes 
the reported incidence Id peak. In this case, for 
which only monthly data are available for both incidence 
and rainfall, not only the rainfall regime is different from 
that in Salvador, but also notice a smaller Pearson's cor- 
relation coefficient (c=0.69) between rainfall and dengue 
incidence (see Figure [6]). 

Mossoro's larger surface of A c = 2110 x 10 6 m [39] di- 
rectly influences spatial parameters, leading to a lattice 
size L = 554. As this incidence data is not georeferenced, 
/i could not be directly estimated. However, taking into 
account that the city is a less developed urban center with 
a smaller population than Salvador (p c =234.390 habi- 
tants (32]), we consider a smaller value of /i = 1.0 x 10~ 6 . 
The values of other parameters were chosen according to 
the already discussed procedures. We observe that the 
normalization factors for the actual (71) and the simu- 
lated (72) incidence time series are very similar. The 
results in Figures [4] and [6] show that the model is robust 
enough to simulate dengue incidence for cities with high 
and low rain intensities, and different Pearson correlation 
coefficients. Thus, such results indicates that, besides 
the importance of periodic forcing, the epidemic behav- 
ior of vector transmitted diseases are heavily dependent 
on entomological and epidemiological aspects that are 
also caught by the model. 

To better understand the forcing effect, the behavior of 
CA model has been followed for large time intervals. We 
consider that the exactly the same rainfall incidence ob- 
tained from one-year pluviometric data is repeated peri- 
odically [12]. Our results indicate the that periodic forc- 
ing leads to modulated responses. However, if we dis- 
allow the possibility of new exogenous infected sources 
(due, e.g., to an infected visitor), the amplitude of the 
epidemic outbursts does not remain the same. If the 
same parameter values as in Figure [4] are used, the re- 
sults in Figure [7] indicate that IHn oscillation amplitude 
reaches its maximum value in the second year, when it 
starts decreasing in a steady way. It is interesting to note 
that, at the same time, the M and A populations do not 
decrease in a similar way. This indicates that, in a closed 
environment, the number of individuals carrying active 
virus and a relatively weak screening effect due to a small 
RH population, turns it difficult to trigger new epidemic 
events. Note that, after five years, the number of suscep- 
tible individuals SH in the population is still very high: 
97% for the parameter set that causes the incidence go 
to zero. 

On the other hand, Figure [7] also shows that changes 
in the parameter values, favoring virus permanence in 
M and H levels for a longer time, may lead to the op- 
posite landscape, with a long period during which the 
yearly amplitude of IHn population increases monoton- 
ically. In such cases, the amplitude decreases only when 
a large fraction of the H population has become infected 
and switched to the RH state. Note that this is not yet 
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FIG. 7: Predicted average IHn for a large time interval of 
5 years as function of entomological features, a) Different 
values of probability of mosquito death: black squares and 
dark-grey triangles indicate, respectively, p m d = 1/7 = 0.143 
and Pmd = 1/8 = 0.125, while the corresponding normaliza- 
tion factors are 15 and 202. b) Different values of human 
viremia period: black squares and dark-grey triangles indi- 
cate, respectively, t v — 5 and t v = 7, with normalization 
factors 15 and 219. Other parameter values are the same as 
in Figure [4] Averages and respective error bars (grey) have 
been taken over M samp i es = 20. 



the situation, after 5 years evolution period, for such al- 
ternative time evolution scenarios. There we still find a 
large fraction of SH susceptible individuals: 52% (see 
Figure [7^i) , where we introduce a variation of probabil- 
ity of mosquito death that is the inverse of expected life 
time of mosquito (p m d = 1/8 = 0.125), and 46% (see Fig- 
ure [7Jd) where the variation occurs in the human viremia 
(t v = 7). The different values of normalized factors in 
both cases indicate how these parameters increases the 
number of IH^. 

This dramatic dependence of the size of successive epi- 
demic events in isolated environments turns to to be a 
unexpected result of our model. As far as we know, this 
effect, resulting from a local interaction between the three 
CA levels, has not been previously discussed in the liter- 
ature. 
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B. The mobility effects: spatiotemporal patterns 

and vector control 3 J 



Spatiotemporal patterns resulting from geo-referenced 
data of the actuald epidemics of Salvador in 1995 have 
been reported elsewhere [35 . They can be compared to 
the CA simulated spatiotemporal patterns, which have 
been generated with the help of the G2 graphic pack- 
age [37] ■ To this purpose, it is necessary to assume that, 
in each CA level, more than one individual can live in 
each lattice site. We consider that the total population 
of the city is represented by the CA cells, assuming the 
inhabitants are a gaussian distributed among the cells 
with a mean value of 50 humans per cell. With this as- 
sumption, the model is able to reproduce qualitatively 
the main features observed in actual spatiotemporal epi- 
demics patterns [55] , 

In Figure [8] we illustrate spatiotemporal patterns for 
A, M and H populations in characteristic time steps. 
For the sake of a better visualization, we choose a small 
value of lattice size (L = 79). As initial condition, we 
assume an infection seed, represented by one IH site in 
the H level. Further, due to a previous large rainfall 
event, the E and SM states of the A and M levels are 
largely populated. From this time on, epidemics starts 
around the site where the seed was located. SM changes 
into EM state, disseminating the disease into other H 
sites, while increasing the radius of the primary epicenter. 
Due to H and M mobility, some secondary epicenters are 
formed. In this case, without any control strategy, the 
epidemics evolves naturally until its end. Figure |8]reveals 
qualitative similarities to the main features presented in 
[35] : the persistence of the epicenter of the epidemics, 
the emergence of secondary epicenters, and an irregular 
shape of each epicenter. 

Secondary epicenters at large distances from the origi- 
nal seed are a direct consequence of the mobility effects, 
which are well accepted to be an important feature for 
dengue transmission urban centers. Indeed, if /i = 0, the 
shown spatiotemporal pattern is replaced by a diffusion- 
like pattern with a single epicenter. However, fi also plays 
an important in reducing time series fluctuations, an ex- 
pected 'mean-field' effect related to the global infection 
probability. This effect is made clear in Figure [9] The 
curves also show that non-zero values of /i introduce a 
time delay effect extending the duration and the inten- 
sity of the epidemics process. Indeed, the large difference 
in the normalization factors for both curves indicates that 
\x is directly related to a much faster epidemic dissemi- 
nation. 

Until today, no efficient vaccine against dengue could 
be devised. Therefore, actions towards vector control 
constitute the only public health policy to reduce the 
deleterious effect of the disease. Even so, there are still 
controversies regarding whether vector control actions 
are more reliable in the A or M phases. As the CA 
model is able to successfully reproduce epidemics data 
and follow the dynamics of the disease for longer periods 
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FIG. 8: Simulated spatiotemporal configurations of cumu- 
lated cases, consider one sample and the following parameter 
values: L = 79 , p m d = 0.143 p em h = Pehm = 0.75 fj, = 0.001, 
t e = 5, ti = 7 , tp — 3 , t; m = 7 , t lh — 5 , t v — 5. Four snap 
shots for each lattice: a) Mosquitoes (M); b) Humans (H); 
and c) Aquatic phase (A). For on line version: ((M): empty 
site - blue, SM - green, EM - grey, IM - red), ((H): SH - 
green, EH - grey, IH - red, RH - blue), ((A): empty site - 
blue, E - green, L/P - grey, B - red). For printed version: 
((M): empty site - white, SM - light-grey, EM - dark-grey, 
IM - black), ((H): SH - white, EH - light-grey, IH - dark- 
grey, RH - black), ((A): empty site - white, E - light-grey, 
L/P - dark-grey, B - black) 



of time, it can also provide useful insights regarding the 
effect produced by different vector control mechanisms. 

To this purpose, let us consider the dependence be- 
tween the epidemic threshold and the vector control pa- 
rameter Pamd- We have performed a large number of 
independent simulations for different values of p am d- We 
evaluated Ith with the help of equation where the 
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FIG. 9: The mobility parameter effect: normalized IHn 
assuming M samp i es = 200 and the same parameters values of 
figure [4] except p that is assumed the following values: 0.0 
(black square) and 0.02 (white circle). The normalization 
factors are, respectively, 17 and 1867. 



time average was replaced by sample averages. Thus, Ith 
is directly related to the probability that an individual 
living the the simulated urban center gets infected within 
a one- year time span. 

The results in Figure 10 show that, when \i = 0, the 
dependence between I t h and p am d follow a power law 
behavior, I t h — ap" md , with large values. It clearly shows 



that effective policies aiming at a reduction of the vector 
reproduction in its own environment produce substantial 
reduction of affected population. This effect is still more 
expressive and relevant when we consider more realistic 
situations, in which human and vector population move 
in the urban space. Indeed, when fi > 0, Ith decays with 
respect to p am d in a faster way the points of fit quite 
well to a stretched exponential I t h = b\ exp[— ^(pfmd)]- 
Moreover, as expected, the epidemic threshold is larger, 
for any value of p am d, when /1 > than when ji = 0. 



V. CONCLUDING REMARKS AND 
PERSPECTIVES 

The three level CA model investigated in this work 
presents several features that allow for a quantitative 
reproduction of actual time series of dengue epidemics. 
Besides the usual local interaction steps based on SEIR 
compartment models, the most important novelties are: 
i) the use of the climatical data as input data; ii) the 
A — M and M — H inter-level interactions; iii) the inclu- 
sion of short-range vector mobility and long-range human 
mobility. 

The model is robust with respect to the range of pa- 
rameters considered in the literature, and to its ability 
in reproducing time series of dengue epidemics in differ- 
ent urban centers. The climatic input data as well as 
the procedure used for estimating the parameter values 



FIG. 10: The vector control analysis: the simulated Ith x 
vector control parameter (p a md)- We consider M samp i es = 20 
and the same parameters values of Figure |4j except for the 
parameter fi — and /1 = 10 -5 . The corresponding values are 
indicated by, respectively, white triangles and black squares. 
The parameter of the power law fitting (grey), for fj, = 0, 
is the exponent a = -2.13 ± 0.04 and a = -1.72 ± 0.05. 
The parameters of the stretched exponential fitting (grey), 
for n = 10~ 5 , are 61 = 1218.04 ± 7.00, b 2 = 86.44 ± 4.00, and 
the exponent (3 = 1.14 ± 0.01. 



are able to catch the diversity of the time series dengue 
incidence for different cities. Although we have mainly 
focused our analysis on the human population, the CA 
model also provides useful insights on the behavior of 
the vector population, which will be presented in a fu- 
ture work. 

The effect of periodic forcing allows us to suggest effec- 
tive measures to reduce the probability of recurrent out- 
bursts. Indeed, the effect of an increased infected vector 
life time is found to be very important to alter of the 
magnitude of epidemic events. 

The analysis of vector control shows that, as expected, 
it indeed produces a decrease in the probability of hu- 
man infection. However, we have shown that this effect 
is more relevant when vector and human mobility are 
taken into account. In this case, the infection probabil- 
ity decreases according to an stretched exponential, while 
a power law behavior is observed when the no mobility 
assumption is taken into account. 

Perspectives for further work on this model are of two 
kind. The first one amounts to investigate the the im- 
pact of different strategies of vector control on dengue 
transmission as well as to discuss the detailed behavior 
of M and A populations subject to those strategies. A 
more ambitions goal is to achieve the quantitative repro- 
duction of spatial patterns. This requires a more pre- 
cise local characterization of spatial units, as well as a 
more precise GPS georeferencing data. This way, the CA 
model can help to plan improved vector control policies 
from the spatial point of view, attacking mainly the most 
important focus for the propagation of the epidemics. 
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Parameter 


Range of values 


Egg period (t e ) J27J 


4-5 days 


Larvae phase period (ti) 27 


5-7 days 


Pupae phase period (t p ) 27 


2-3 days 


Latent period of virus in the mosquito (t; m ) 9, 22 27, 29 


7-20 days 


Latent period of virus in the human (tih) 9, 22, 29 


2-12 days 


Viremia period (t v ) |9| 1221 1301 


3-7 days 


Probability of transmission human-mosquito (p e hm) 1311 


0.5-1.0 


Probability of transmission mosquito- human (p ern h) [H] 


0.5-1.0 


Probability of mosquito death (p m d) [9, 25, 32 


0.128-0.25 



TABLE I: The parameter range of values of temporal parameters and the probabilities of transmission H-M and M-H, and of 
death mosquito according to the literature. The baseline values were chosen for the simulations of the model. 



